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THE INFORMATION IN THIS SOFTWARE IS SUBJECT TO CHANGE WITHOUT NOTICE 
AND SHOULD NOT BE CONSTRUED AS A COMMITMENT BY DIGITAL EQUIPMENT 
CORPORATION. 

DIGITAL ASSUMES NO tery: 


IBILITY FOR THE USE OR RELIABILITY OF ITS 
SOFTWARE ON EQUIPMENT WHICH IS NOT SUPPLIED BY DIGITAL. 
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; facility: 
VAX/VMS PL1 runtime Library. 
: ; abstract: 


Runtime routine performs fixed decimal (packed decimal) division. 
The routine is called when precision and scale Pogutrenents for 
the quotient imply age precision division. 

only called when such mu 

when the divisor has a precision of less 


he routine is 
tiple precision divistgn is required and 
0 decimal _tteite. 
(Call pli$div_pk_long if multiple preg sien division 
0 or 31 decimal digits). 


PLISDIV_PKSHORT 


required and the divisor has precision 
; author: Peter Baum 20-jun-1980 
; modifications: 


COOoOOCCOCCOOSOSCSOOOOSOOOCSOSOSOOSOOCOOSCOOOOOOOOOO 
SOOSCSCOCSOOSOSOOCOOOOOOOOOOSOOOOOOOOOOOSOOOOOCOOSOOoOO 


1-002 Bill Matthews 2?9-September-1982 
Invoke macros $defdat and rtshare instead of $defopr and share. 
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: routine: 
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PLISDIV_PKSHORT 16*SEP=1984 02:22:40 VAX/VMS Macro v04-00 Pa | 
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documentation file: Cpl1.doc.codegen] THEORY.MEM 
functional description: 
This routine calculates: 


z=ex/y | 
let a = scale(z) + scale(y) - scale(x) - 31 + ereexm 
b = scale(z) + scale(y) = scale(x) + prec(x 
c = 31 = prec(x) 
d = 31 = prec(y) 


this routine is called if b > 31 and d > 1 


Prior to the call: 
if c not 0 then shift x left by c. 
Thus x is a 31 digit packed decimal. 


input: 


address of dividend (shifted left by c) 

address of divisor 
) precision of divisor (high order bytes zeroed) 
) address of quotient 
) precision of quotient (high order bytes zeroed) 
) a as defined above(high order bytes zeroed) 
) das defined above(high order bytes zeroed) 


MmMPnN——@LfoO 


output: . eS 
quotient returned at address specified by 16(ap) 


| 

| 

| 

# of arguments 
variable usage: | 


size 
in 
variable digits use 


x (ap) 31 Dividend 
y\ep), py(ap) Divisor 
a 


SOOSCOCSCOSOOOOOOOOOOSOOOOSOOSOOOSOOSOOOSOOOOOCOOOOOOOOOOOOOOOOOOOOOOOo 
SOOOCSCSOSOSOSCSSOOOSOSOSOOCSOOSOOSOOSCOOSOOSOSOOCSOSOSOSOSOOSOSOSOOSOOSO OOOO OOOOOSOOOOO 
COOCCSOCOOSOOOCOSOSOOOCOCOOOO SOOO SOSOOCOCSOOOOSCOSOSOCOSOOSOOSOOOOOOOOCOOOOSOOOOOOoOO 


Dm Tm Je |e a a | we | ee fe ee) we eee) we) we we ww ow wl w lolol olololeloleoleleleolelelolelelelelolelelolelelelelolol ola a) 
= = |§ DOO DOGD9900000O0 OOO OOO 09 00 09 09 09 09 09 09 09 09 SIN NN NINN NIN NNO AAAAAAAAOOUUN 
FWD ODO OONAU EWN 9 OONAU EWN 9 OD NAME WN 0 OONOAU SW OO ODNOULS WN —O0OO 


SO a as a a 2) — 4 2 2 


by p ecccee Binary number that gives precision of y 

z(ap pz(ap) Quotient . 

pz(ap) ------ Binary number that gives precision of 2 

(sp) 31 Initially abs(x); successive remainders as 
algorithm progresses. 

stkz2(sp) d Temporarily holds the next d 
digits of quotient. : 

stktl(sp) 31 Temporary because packed instructions 


PLISDIV_PKSHORT 
1-002 


=" 
an 
te 
nom 
mm 
vv 
te 
—— 
ono 


Bo 99:85:10 EPCINTL. SREIPLIPROIVS.MAR; 1 


don't allow overlapped operands 
stky(sp) 31 Holds abs(y) 
stksign(sp) 2 bits used to indicate the sign of the 
quotient. 00=+, 10=+, Ol=-; via incb 


register usage: 
register use 


r a = additional digits of precision required beyond prec(x) 
r stky(sp) = address of divisor 
r8 py(ap) = precision of divisor 
r9 r = number of additional digits of the quotient 
10 i to be found for next step 
r z(ap 
ri d = 31 = prec(y) = max. no of digits obtained each iteration 


optimization notes: 


1) Optimized for speed, not space. 

2) Optimized for y > 0. ; 

3) Assumes speed for register to register operations are the same 
for byte operations and longword operations. 

4) Many packed instruction ei were timed. Do not change 
unless actual tests are made to determine relative speed. 
Tests were made on 11/780 and Comet. 


SOOSCOSCOSCOOOOSOOOOOOOOOOOOOOOOSOSOOSOSOOOOOOCOOOOCOOOOOOOCOOOOOOOOOO 
SIDS DS DDS DS DS DDD DS TTT BB ES SB BB AAW nonnononononononong) 2] 2 2 
AO OD NAW EWN HO OOO NAME WN O OO NAU EWN O OD NOAU EWN 0 DONAU EWN" OOONOU 


SOSSCSOSOOCSOSOSSSOOSOSOOSOOCOOOOSOOSOSOSCOOSOOOSoSO 


SOOO SOSCSOOOSOSOOCOOOOOOOSOOSOOSOOoOSoOSoOo 


SOCSCCOSCOOCOOOOCOCOOOCOOOOOOoOoOoOoOoO 
a ee ed dd dd 


stack offsets for work area 
Soffset 0,,<- 


<,16>,- sabs(x), 31 digits 
<stkz2,16>,- :z2 31 digits 
<stkt1,16>,- 3t1 31 digits 

<stky, 16>,- sabs(y) ; ; 
<stksign,{>,- :sign of quotient, 2 bits 
<stklen,O>,- slength of work area 


: parameter offsets 


$Soffset 4,,<- : hae 

<x?,> :x = dividend by reference 
<y>,> zy = divisor by reference 
> ae sprec(y) by value 

<z>,° iz = papt Sens by reference 
<p2z>,- zprec(z) by value 
<consta>,- za by value 

<consté>,© id by value 


PLISDIV_PKSHORT 16-SEP-1984 :40 VAX/VMS o V04-00 Page 4. PL 
1-002 6-SEP=1 198 94: $6 15 LOE TRIL L.S SPLIPKDIVS.MAR; 1 (1). 1- 
000 172 ; 
tn 
rtshare 
0000 175 ; 
0000 178 3 constant data area 
4 7 3 
0c 4! 178 zero: packed +0 slocal packed decimal constant zero 
449 39 ; local symbol definitions 
0000000F i498 186 bytes_to_sign=15 bytes to sign for fixed decimal 31 
$001 184 3 
0001 185 ; 
0001 186 
CFFC 8h Ht entry pli$div_pkshort,*M<iv,dv.r2.r3,r4,r5,r6,r7,°8,r9,r10,r11> 
it 13) sinitialize registers and temporaries 
5E BF AE 9E€ 0003 191 movab -stklen(sp),sp smake room for temporaries 
5A 10 AC 00 0007 136 mov z(ap),r10 isave address of quotient 
57 30 AE 9E€ OO0B 19 movab stky(sp),r7 saddress of divisor 
58 OC AC DO OOOF 194 movl py(ap),r egrecisten of divisor 
5B 1C AC 00 0013 195 movl constd(ap),r11 = 31 - prec(y) 
56 18 AC DO 0017 196 mov consta(ap) ,r6 + = _scale(z) + scale(y) - scale(x) | 
001B 197 3 > 31 + prec(x) 
40 AE 94 QO1B 198 clrb sttsiontse) iclear sign flag 
6— 04 6C 1F 34 QOIE 199 movp #31,ax(ap), (sp) smove x, set cond. code 
26 14 0023 200 botr 50$ sbranch if x>0 
1E 19 0025 201 blss 40$ tbranch if x<0 
0027 sog : 
0027 203 ;x = 0 
0027 she : 
08 BC 58 ODS AF 00 37 0027 05 ns ag #0,zero,r8,ay(ap) spot condition code 
OA 13 OO02E 206 beq branch if divide by 0 
14 AC 00 CB AF 00 ei F8 Oo 3 207 ashp #0,#0,zero,#0,pz(ap), cr16) 32 = 0 | 
04 0039 208 ret 
BF AF OO C2 AF O00 27 QO3A 209 308: divp #0,zero,#0,zero,pz(ap),(r10) ;cause divide by 0 
6A 14 AC 0041 
04 0044 $19 ret 3 | 
0045 11 ; : : 
0045 12 ;x not 0, determine sign of x 
Ope e 13 3 
45 214 os: | 
40 AE 96 0045 215 incb stksign(sp) zset low order bit 
OF AE 97 et gig | decb eng, to *Bign(sp) zx <0 so make it positive | 
0048 18 sdeternine si n 
0048 19 sy may be ebis” point 
0048 0 : code RRA FS for y>0 
0048 1 
67 08 BC 58 34 ints 4 S08: 8 ,ay(ap),(r7) into t 
mov r a r smove nto temporary 
oc 8618 $50 4 bgea 60$ _ sbronch {¢ y>=0 
40 a 3 b2¢ 5 incb stksign(sp) -set neg indicates 
AG AF OO O&8 BC 5 2 05 6 subp6 r8,aylap),#0,zero,r8,(r7) sconvert to positive 


E 10 


‘start of multiple precision divide 


1F 10AE SB 67 458 = 25 100$:  mulp r8,(r7),R11,stkz2(sp) ,#31,stktl(sp) ;tl=y*z2 


PLISDIV_PKSHORT 16-SEP-1984 74 AX/VMS o V04-00 Page 5 
He hy 6-SEP=1 138¢ oF: $6; 18 Ler TRIE L.S SPLIPKD 1VS.MAR;1 . (1) 
67 05D 
OSE 60$: 
Q5E ; : 
bee pstart of divide proper; setup 
=.=. sme ae $08 cmpp4 #31,(sp),r8,(r7) px<y? 
13. «19 One blss 95$ sbranch if x<y, i.e. shift of d is o.k. 
0065 sy <x 
0065 ; 
6A 14 AC E 67 58 7 0065 divp r8,(r7),#31,(sp), ath ap), (r10) 
1F 6A 14 AC 67 58 5 8860 mulp r8,(r7) spz(ap), (r10), #34 siuttinet Ttle(x/y)ey 
20 AE 074 
14 11 0076 brb 110$ 3 
0078 3 
BEB nc 
14 AC 00 83 AF 00 rn F8 444 $5s: ashp #0,#0,zero,40,pz(ap),(r10) ;clear quotient 
11. 11 0081 brb 115$ : 
0083 
0083 
0083 
0083 
20 AE 008A 
6— 1F 20 AE IF $6 008C 110$: subp4 #31,stktl(sp),431,(sp) ;x=x-t1 
41 009¢ beq 150§ tbranch if remainder = 0 
0094 idetermine r, the number of the next low order digits to obtain 
59 5B DO 0094 415$: movil = r11,r9 ir=d 
5B 56 Di 0097 cmpl r6,r11 ya>d? 
03 14 OO9A bgtr 1308 corenah if larger 
59 56 00 009C mov l r9 
20 AE 1F OO 6€ +1F 59 =&§F8 QO9F 130$: ashp 5°41, (sp),#0, “2 sstkt (sp) spurs A bg by r 
6— 20 AE 1F 34 OOA7 movp 8 31 »stktl (sp) opy ba in 
14 AC 00 6A 14 acon + F8 Beat ashp r9 »pz(ap), (F160). #0, eaten). st t1isp) rghit z left by r 
6A 20 AE 14 AC 34 00B6 movp pz (ap) sthst (sp). (r10) scopy back into z 
10 AE SB 6€ IF 67 58 27 O0BC divp §.#31, (sp) rll,stkze(sp) ;22(d)=x/y 
6A 14 AC 10 AE 5B 20 0004 addp4 rif sgtkz2sp), pz(ap).(r10) ;z=z+z2 
56 59 C2 00CB subl2 9,r6 jazacr 
os. Tz bneq 106$ sbranch if more 
16 40 AE fe aioe stksign(sp),155$ tbranch if quotient <0 
re by 


tremainder = 0 


150$: remainder = 0 
ashp r6,pz(ap),(r10) 40, paten!. stktl(sp) ;account for scale 


blbs athe tentap) 160$ sbranch if quotient < 0 
og pz(ap),st ti(sp), (r10) ; Copy back inte quotient 
re $ 


0D 40 A 
6A 20 AE 14 AC 4 


iquotient <0 
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PLISDIV, PKSHORT 


20 AE 6A 14 AC 34 
FFO7 CF «600 20 AE 14 AC 23 
6A 14 AC 


movp 
subp6 
ret 
-end 


"SSEp=198e 9:95:70 EBCIRTC. SRETPL LPKDIVS.maR; 1 


pz(ap),(r10),stktl(sp) 
zen 
pz(ap),stktl(sp),#0,zero, eatae!, Mt fd) 


poner Rs ew: into temp 


holds quotient 
smake z negative 


Page 
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PLISDIV_PKSHORT 
Symbol fable 


6483 byies (1 


= 


BYTES _TO_SIGN = OOOOF 
CONSTA ~ 8 O08 
CONSTD 000001C 
DIR... = 00000001 
PLISDIV_PKSHORT se RG 02 
PY 000000C 
4 
STKLEN 0000041 
STKSIGN 00000040 
STKT1 44 34 
STKY 00000030 
STKZ2 aes ah 
x 00000004 
Y 00000008 
Z 00000010 
ZERO 90000000 k 02 
wee eB eB BeBe we ee eS + 
: Psect synopsis ; 
PSECT name Allocation PSECT No. Attributes 
S 00000000 < 0.) 00 ¢ O.) NOPIC USR CON 
S$ 00000041 (¢ 65.) 01 ¢ 1.) \NOPIC USR_- CON 
PLISCODE OO0O0000FD (¢ 253.) 02 ¢ 2.) PIC USR CON 
(deanna 
! ;_Performance indicators ! 
Phase Page faults CPU Time Elapsed Time 
Initialization 9 00:00:00.07 00:00:00.32 
Command processing 77 00:00:00.51 00:00:01.69 
Pass 67 00:00:01.18 00:00:02.27 
Symbol table sort 0 00:00:00.01 00:00:00.01 
Pass 51 00:00:00.63 00:00:01.16 
Symbol table output 2 00:00:00.02 88:88:00 B5 
Psect synopsis output 1 00:00:00.02 00:00:00.0 
Cross-reference output 0 00:00:00.00 00:00:00.00 
Assembler run totals 207 00:06:02.45 00:00:05.49 


The shereses 2 set Limit was 750 pages. 

pages) of virtual memory were used to buffer the intermediate c 
There were 10 pages of symbol table space allocated to hold 16 non-local and "3° lecat symbols. 
83 source Lines were read in Pass 1, 
pages of virtual memory were used to define 


12 object records in Pass 
macros. 


producin 


AX/VM 
PLIRT 


“1986 19:85:10 


o V04-00 
JPLIP KDIVS.MAR; if 


ete NOEXE NORD _ NOVEC BYTE 
NOSHR EXE RD 


SHR 


EXE 


WRT NOVEC BYTE 
RD NOURT NOVEC LONG 


1-| 


| 
| 


H 10 
VAKetd Recto Run Statistics "B28 : Pa] 984 4: 96: 8 Oe tR TL E58. Yer DIVS.MAR; 1 iva Gy 
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Macro Library name Macros defined 
3 $2558DUA28: PLIRTL .08 JPLIRTMAC.MLB; 1 3 
$255$DUA SYSLIBI STAR LET.MLB;2 ° 
TOTALS (all Libraries) 


44 GETS were required to define 3 macros. 


There were no errors, warnings or information messages. 
MACRO/ENABLE=SUPPRESSION/DISABLE=TRACEBACK/LIS=LIS$:PLIPKDIVS/OBJ=OBJ$:PLIPKDIVS MSRC$:PLIPKDIVS/UPDATE=(ENH$:PLIPKDIVS)+LIB$:PLIRTM 
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